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We analyze a scheme to manipulate quantum states of neutral atoms at individual sites of op- 
tical lattices using focused laser beams. Spatial distributions of focused laser intensities induce 
position-dependent energy shifts of hyperfine states, which, combined with microwave radiation, 
allow selective manipulation of quantum states of individual target atoms. We show that various 
errors in the manipulation process are suppressed below 10 -4 with properly chosen microwave pulse 
sequences and laser parameters. A similar idea is also applied to measure quantum states of single 
atoms in optical lattices. 

PACS numbers: 03.67.Lx, 39.25.+k, 03.75.Lm,03.75.Mn, 



I. INTRODUCTION 

Neutral atoms trapped in optical lattices are excellent 
candidates for quantum computation because they are 
well isolated from environment, leading to long coherence 
times, and easy to use for storing and processing quantum 
information (HOI • In optical lattices, controlled interac- 
tions between atoms may be implemented effectively, and 
highly entangled states of many atoms may be created in 
a single operational step Uil- Lil- LJ| • More importantly, 
very efficient schemes for quantum error corrections |(j 
and fault-tolerant computing can be straightforwardly 
implemented because of parallel operations in optical lat- 
tices [S| and long coherence times. 

An outstanding challenge for quantum computation in 
optical lattices is the selective manipulation and measure- 
ment of quantum states of single atoms, because spatial 
periods of typical optical lattices are shorter than op- 
tical resolutions (long wavelength lattices may provide 
single site addressability, but it is difficult to implement 
Mott insulating states with one atom per site necessary 
for quantum computation in such lattices Single 
qubit operation is not only a building block of universal 
quantum computation |10|. but also an essential ingre- 
dient for the recently proposed one way quantum com- 
putation, where quantum information is processed by 
performing single qubit rotations and measurements on 
entangled cluster states 0. Experimentally, entangled 
cluster states have been realized for neutral atoms by 
using controlled cold collisions in spin-dependent optical 
lattices Therefore, implementations of single atom 
manipulation and measurement may eventually lead to 
universal quantum computation in optical lattices. 

In this paper, we show that high fidelity selective ma- 
nipulation and measurement of single atoms in optical 
lattices can be achieved with the assistance of focused 
lasers 12]. Consider a deep two dimensional optical lat- 
tice with one atom per lattice site . The logical qubit 
basis of each atom is formed by two hyperfine states 
that can be coupled coherently using microwave radia- 
tion. The couplings are same for all atoms because of 



the degeneracy of their hyperfine splittings. In the pres- 
ence of spatially varying external fields, the degeneracy 
may be lifted and atoms separated by certain distance 
can be individually manipulated. For instance, individ- 
ual atoms separated around 2.5/im have been selectively 
addressed in an experiment using magnetic field gradi- 
ents , but such a method is not useful for atoms sep- 
arated by a half-wavelength of typical optical lattices, 
where impractically large gradients or small fluctuations 
of the magnetic fields would be required. 

The degeneracy of hyperfine splittings may also be 
lifted in the presence of focused lasers that induce 
position-dependent energy shifts of hyperfine states (light 
shifts) [l5| due to the spatial distributions of laser in- 
tensities. Through varying intensities and detunings of 
the focused lasers, we may adjust the difference of hy- 
perfine splittings between neighboring atoms to a regime 
where microwave radiation affects mainly target atoms, 
while impact on non-target atoms is strongly suppressed 
with properly chosen microwave pulse sequences. Vari- 
ous errors in the manipulation process are found to be 
below 10 except for the spontaneous emission proba- 
bility of atoms in the focused laser. Finally, we show that 
the quantum states of individual atoms in optical lattices 
may also be measured using position-dependent energy 
shifts of hyperfine states. 

II. SELECTIVE MANIPULATION OF SINGLE 
ATOMS 

We develop the scheme as realistically as possible, and 
illustrate it using 87 Rb atoms, although the technique is 
applicable to other species as well. Consider 87 Rb atoms 
confined in a two dimensional (xy plane) optical lattice 
with wavelength A = 850nm. Perpendicular to the plane, 
the atomic dynamics are frozen out by high frequency op- 
tical traps [16j . The lattice is ramped up adiabatically to 
a potential depth Vl — 5QE r such that the Bose-Einstein 
condensate is converted into a Mott insulating state with 
one atom per lattice site p|. Here E r denotes the re- 
coil energy E r = h 2 k 2 /2m « h x 2tt x i.lSKHz. Two 
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Figure 1: (color online). The scheme for manipulating single 
atoms in optical lattices using focused lasers and microwave 
radiation, (a) S7 Rb hyperfine structure. Ao is the detuning 
of the focused laser from the transition 5 2 S 1 / 2 — > 6 2 P 3 / 2 - (b) 
Position-dependent hyperfine splitting induced by the spatial 
distribution of the focused laser intensity. Microwave radia- 
tion couples two hyperfine states of atoms, (c) Optical poten- 
tials for atoms around the focused laser. The length unit is 
the wavelength A = 850nm of the optical lattice. The energy 
unit is the recoil energy E r . Top: a 2D plot for atoms at state 
|0). Bottom: a ID plot along y = 0. Solid and dashed lines 
for states |0) and |1) respectively. 



hyperfine ground states \F = l,m,F = —l) = |0), and 
\F = 2,mF = — 2} = |1) are chosen as the logical basis 
of a single atom qubit at each site. 



A. Position-dependent hyperfine splittings 

To manipulate a target atom A, we adiabatically turn 
on a focused laser that propagates along z axis having the 
maximal intensity located at A (Fig.l). The focused laser 
is obtained by passing an initial large Gaussian beam 
with waist w through a lens with diameter D and focal 
length /, and its intensity shows an Airy pattern I (r) = 

W (.// r ' J o (kfr'r/y/r 2 + / 2 ) cxp (-r' 2 /w 2 ) drf , 
where kf — 27r/A/ is the wavevector of the focused 
laser, Jo i r ) is the zero order Bessel function, and 



G = J Q /z r'J (0)exp(-r' 2 /w 2 )dr'. With properly 
chosen lens parameters (see Tab I), the initial Gaussian 
beam is focused to the diffraction limit, where the 
intensity pattern is accurately described by above Airy 
function. In contrast, commonly used Gaussian beam 
approximation is not precise in such limit. In addition, 
the Airy pattern has a narrower waist for the central 
intensity distribution, leading to a smaller focused laser 
intensity on neighboring atoms, although the pattern 
extends more broadly than a Gaussian distribution that 
has an exponential decay. 



AE i (r) = ^I(r) 



(1) 



for hyperfine ground states |0) and |1), where c is the 
speed of light, Tj is the decay rate of excited state \j), 
Cij is the transition coefficient, u>ij is the frequency, and 
Ay is the detuning of the focused laser for the transi- 
tion |z) — * \j). Different polarizations and detunings of 
the focused laser yield different shifts of hyperfine split- 
tings \AE (r)| = |A£a (r) - AE (r)| between two qubit 
states. Here we choose a <r + -polarized laser that drives 
the 5S — ► 6P transition 18] to obtain a small diffraction 
limit (Fig. 1(a)). The laser induces a red-detuned trap for 
state |0), but a blue-detuned trap for state |1). The wave- 
length A/ « 421nm (i.e. the detuning Ao) is optimized 
to obtain the maximal ratio between energy splittingof 
two qubit states and the spontaneous scattering rate [j] . 

Because of the inhomogeneity of the focused laser in- 
tensity, |Ai?(r)| reaches a maximum at the target atom 
A and decreases dramatically at neighboring sites, as 
shown in Fig. 1. Therefore the degeneracy of hyperfine 
splittings between different atoms is lifted (Fig. 1(b)). In 
Fig. 1(c), the optical potential V for atoms around the 
focused laser is plotted. We see that the minimal poten- 
tial barrier, AV^ 20Er, occurs for atom B at state |0). 
The tunneling rate of atom B to the neighboring site A 
is vj w 2irJ 2 /HE g using Fermi's Golden rule, where J 
is the hopping matrix element and E g is the energy gap 
between the initial and final states. In a symmetric opti- 
cal lattice, Eg is just the on-site interaction between two 
atoms. The asymmetry of sites A and B yields an energy 
gap Eg on the order of the trapping frequency of site A, 
which strongly suppresses the tunneling rate. A rough 
estimate shows that the tunneling time 1/w is about 13s 
for the parameters in Tab. I. 

To avoid excitations of atoms to higher bands of opti- 
cal lattices, the rising speed of the focused laser intensity 
should satisfy the adiabatic condition H\dAEi (0) /dt\ = 
^\co eg \ 2 /\{^> e \dV/dAE l (0)\<Pg)\ 1 where uo eg is the en- 
ergy gap, and |$ e ) are the wavefunctions of the 
ground and excited states, and £ <C 1 is a parameter 



Table I: Experimental parameters. 5 (A/2) is the difference of 
the hyperfine splittings between atoms A and B. For other 
atoms, 8 is slightly larger, r is the spontaneous emission 
probability for atoms in the focused laser during the single 
qubit manipulation process. Pf is the power of the focused 
laser. The recoil energy E r = h x 2n x 3.18KHz. 
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that determines the degree of adiabaticity. Because of 
the high potential barrier for each atom, the wavefunc- 
tions and energy gap may be obtained using the harmonic 
oscillator approximation, and we find 



dAEj(0)/E r 
dt 



^^ritpl-^^)^ Where °0 = V h2 / mE r, ^ M IS 

the potential barrier for atom at position r with state 
and we have approximated the Airy pattern of the 
focused laser intensity using a Gaussian function with 
an effective beam waist w. The rising speed is lim- 
ited by atom B at state |0) that has minimal poten- 

3.6 x 



rfA£ o (0)/g r 
dt 



tial barrier, and is estimated to be 

10 6 £ (AV (A/2) /£,.) 5/4 Hz. The total ramping up time is 
found to be 57/xs for £ = 0.005, which corresponds to a 
10~ 4 probability for excitation to higher bands. 



B. Single qubit rotation 

The position-dependent hyperfine splittings induced 
by the focused laser, combined with microwave radiation, 
can be used to perform arbitrary single qubit unitary op- 
erations on the target atom A. The microwave frequency 
is chosen to be resonant with the hyperfine splitting of 
A, and has a S (r) = (\AE (0)\ - \AE (r)\) /% detuning 
for non-target atoms (Fig. 1(b)). The coupling between 
two states |0) and |1) is described by the Rabi equation 




e- i m{t)/2 
e' x Q(t)/2 -5 




(2) 



where fl is the Rabi frequency and x is the phase of 
the microwave. The evolution of the quantum states is 
equivalent to the rotation of a spin with components S z = 
\ci\ — | Co I , S x = 2 cos£? |ciCo|, S y = 2sin6>|cico on a 
Bloch sphere ^5|, where 9 is the relative phase between 
two states. 5 is the rotation frequency along S z axis, 
while f2 (t) is the frequency along an axis on xy plane 
whose direction is determined by phase For instance, 
X = and tt/2 correspond to rotations along axes S x 
and S y respectively, and the combination of them may 
implement arbitrary single qubit unitary operations 
In this paper we focus on x = although similar results 
for x = tt/2 may be obtained straightforwardly. 

To achieve a reliable single qubit rotation with negli- 
gible impact on non-target atoms, the microwave Rabi 
pulses need to satisfy the following criteria: (i) the vari- 
ations of the occupation probabilities of state 1 1} of non- 
target atoms are extremely small (below 10~ 4 ); (ii) the 
relative phase 9 between two qubit states must return to 
the initial value after the microwave pulses. In the fol- 
lowing, we show how a single qubit rotation on the tar- 
get atom A satisfying these criteria is achieved by using 
Gaussian shaped rotation pulses together with refocusing 
pulses. 

A Gaussian shaped pulse (t) — £1 exp (— Wq£ 2 ) 
(— tf < t < tf) is used to perform a single qubit rota- 
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Figure 2: (color online), (a) Time evolution of A|ci(t)| 2 
during the microwave pulse, fio/i^o ~ l-^l is the Rabi 
amplitude for a pulse whose area is slightly larger than 7r. 
Atoms A (dashed dotted) and B (solid) with initial state 

(co,ci) = a toms A (dashed) and B (dotted) 

with initial state (1,0). (b) A|ci (t/) | 2 is the population varia- 
tion of atom A after the microwave pulse. Initial states (1, 0) 

(dashed dotted); (\fh,i\fh) (dotted), (,/ / |, J\ + i*f± 



(dashed), and 



(solid). 



tion on the target atom A. In the frequency domain, 
the Fourier transformation Q (lo) of the Rabi frequency 
SI (t) of this type of pulse shows a Gaussian shaped de- 
cay with respect to the detuning 5 from the microwave 
frequency. Because of their large detuning S, non-target 
atoms undergo Rabi oscillations with small frequencies 
O (5), which strongly suppress the variations of their oc- 
cupation probabilities. 

This scenario is confirmed by numerically integrat- 
ing Eq. J5J) with different initial states and calculat- 
ing the variation A|ci(i)| 2 = |ci(t)| 2 — l c i(^/)| 2 of the 
occupation probabilities at state |1). A small constant 
ujo = (-|) for the microwave is chosen to avoid large 
impact on non-target atoms, and different pulse areas 
are implemented by varying the pulse amplitude f2o> in- 
stead of the pulse period 2tf. In Fig. 2(a), A|ci(i)| 2 
of atoms A and B are plotted with respect to the scaled 
time u>ot. We see large variations of |ci | for atom A after 
the pulse, which correspond to a single qubit rotation. In 
comparison, there are no obvious changes of |ci| for the 
neighboring atom B. In Fig. 2(b), the total changes of 
\c\\" for the target atom A at time tf are plotted with 
respect to Sl / ■ We see different pulse areas can be 
achieved by adjusting f^o- 

For the target atom (S — 0), the time evolution of 
|ci (t)\ can be obtained analytically by solving Eq. J2J, 
which yields 



|ci (t)\ 2 = (1 - psin fafio/wo + 4>)) /2. 



(3) 



Here p = yf 1 - 4 |c {-t f )\ 2 | Cl {-t f )\ 2 cos 2 9 (-t f ), 
i] = ujq f_ t exp (— cjq<' 2 ) dt' and 4> = 
arcsin ^1 — 2 |ci (— tf)\ 2 ^j /pj . This expression is 
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Figure 3: Manipulation error e = |jci (t/)| 2 — |ci (— t/)| I for 
the nearest neighboring atom _B. e is smaller for other atoms 
because they have larger detuning 5 (r). —7 < coot < 7 in the 
numerical integration of Eq. Initial states are same as 

those in Fig. 2(b). 



in very good agreement with the numerical results 
presented in Fig. 2. In Fig. 3, we plot the single qubit 
manipulation errors for the nearest neighboring atom 
B with respect to VLq/loq. The errors are smaller than 
1CP 4 for different Rabi pulse areas and different initial 
states. We have also performed numerical simulations 
for many other initial states, finding that errors below 
1CP 4 are always achieved. The manipulation errors are 
even smaller for all other atoms because they have larger 
detuning S (r) . 

The relative phase 9 is harder to control than the oc- 
cupation probability because it can not only vary from 
to 27T many times with large frequency 8 (r) during 
the Rabi pulses, but also be easily affected by interac- 
tions with the environment that lead to dephasing. To 
eliminate the variation of 9, we use the following pulse 
sequence that is similar to the refocusing process in NMR 
studies To realize an area a pulse, we adiabatically 
ramp up the focused laser, apply a a/2 pulse on the tar- 
get atom A (it induces a rotation around S x axis) using 
the Gaussian shaped pulse described above, then adia- 
batically ramp down the focused laser in step (i). In this 
step, atoms obtain phase variations A9 determined by 
their detunings S (r) and the dephasing process. In step 
(ii), we apply a fast resonant 7r pulse to all atoms that 
corresponds to an angle it rotation around S x axis. This 
pulse is called the refocusing pulse. In step (iii) we repeat 
step (i), and in step (iv) we apply another tt refocusing 
pulse. The two refocusing pulses do not affect the pulse 
area for target atoms because the combination of them 
corresponds to a 2-7T rotation around S x axis. However, 
the effective directions of the rotations around S z axis for 
all atoms in step (iii) are reversed from that in step (i) 
by two refocusing pulses, therefore the phase variation 
in step (iii) becomes — A9. These two phase variations 
cancel each other and the relative phase 9 returns to its 
initial value after the pulse sequence. 



C. Fidelity analysis 

In an experiment, a difficult parameter to specify is the 
spatial distribution of the focused laser intensity that de- 
termines the detuning 5 (r) , the central parameter of the 
scheme. However, our scheme only requires a large S (r), 
not any specific value, therefore this is not a significant 
difficulty. Another important issue in the experiment is 
the misalignment of the focused laser from the minimum 
of the optical lattice potential. A small displacement Ace 
of the focused laser induces a detuning of the microwave 
scaling as (Ax) 2 from the hyperfine splitting between two 
qubit states of the target atom, and thus reduces the fi- 
delity of the single qubit rotation. For Ax — lnm, we 
estimate the detuning to be 2-7T x 2>Hz and the fidelity of 
the single qubit rotation is degraded by 2 x 10~ 4 . Gen- 
erally, errors due to small mis-detuning of the microwave 
may be corrected using composite pulses technology de- 
veloped in the NMR quantum computation [Tgj. 

In the single qubit manipulation, the probability of 
spontaneous scattering one photon for atoms in the fo- 
cused laser is estimated to be r — 6 x 10~ 4 , which is the 
worst parameter in the scheme. This parameter is lim- 
ited by the need for a maximal ratio between the vector 
light shifts and the spontaneous scattering rate, which 
does not allow arbitrarily large detunings. 



III. SELECTIVE DETECTION OF SINGLE 
ATOMS 

The energy shift induced by the focused laser is around 
several hundreds KHz, which is much smaller than the 
decay rate T of excited states and cannot be used to selec- 
tively measure quantum states of target atoms. To obtain 
large energy shifts, we transfer the target atom to other 
magnetic sublevels (|0) -> |1) = \F = 2,m F = 2), |1) -> 
|0)) using several tt microwave pulses with the assistance 
of focused lasers. Then a focused cr + -polarized laser reso- 
nant with the transition |1) — * |6 2 P 3 /2 : F' — 3, m! F , = 3) 
is applied to detect atoms at state |1). In this process, 
resonant fluorescence is observed if and only if the ini- 
tial state of the target atom is |0) and the probabil- 
ity to detect an atom yields |co| of the target atom. 
The detection laser can be focused on the target atoms 
as well as the focused laser beam using for single atom 
manipulation because they both use 55 — » 6P transi- 
tion. In the experiment, a constant 30G magnetic field 
along x-axis may be used to induce about 2-7T x 84MHz 
Zeeman splitting between states |1) and 1), while 
2tt x 111.6MHz between |6 2 P 3/2 : F' = 3, m' F , = -l) 



and 6 2 P 3 /2 : F' = 3,m' F 



3/2 

tection laser is about 
the transition |1) — ► |6 2 P- 



= 3) [l5j ■ Therefore the de- 
2tt x 27. 6M Hz detuned from 
1} -> |6 2 P 3/2 : F' = 3,m' F , = -l), which 
yields a maximal ratio (T 1 /2S 1 ) 2 I (r) // (0) < 2 x 1CT 5 
|15| between the photons scattering on non-target and 
target atoms, that is, 5 x 10 4 photons have been scattered 
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on the target atoms to induce one photon scattering on 
the neighboring non-target atoms. The impact on non- 
target atoms in the detection process may therefore be 
neglected. 

IV. DISCUSSION AND CONCLUSION 

We emphasize that such high fidelity single atom ma- 
nipulation cannot be accomplished using magnetic field 
gradients. To obtain the same detuning 6 (r) as that 
using focused lasers, an impractical magnetic field gra- 
dient (~ 3.6T/cm) is required. While in a typical gra- 
dient lOG/cm, a single atom rotation takes about 40ms 
(much longer than the cold collision time (~ 100/zs) for 
many-qubit gates 0,0), and requires unattainably small 
magnetic field fluctuations (~ fO~ 6 G). In addition, the 
magnetic field gradient method can only be used to se- 
lectively address atoms along one dimension, while the 
focused laser scheme can select two dimensions. 

We notice that single qubit rotations may also be per- 
formed using two-photon Raman transitions [20| with the 
assistance of focused lasers. In this case, the beam waists 
of the Raman pulses need to be relatively large (~ 2/im) 
so that small misalignments of the lasers do not cause 
large changes of the laser intensities on the target atoms 
that may diminish the fidelity of single qubit rotations. 
Because two-photon Raman transitions only affect atoms 



inside the Rabi pulses, atoms separated by a long dis- 
tance may be addressed simultaneously with different 
pulse sequences, which may reduce the total computa- 
tion time. 

In summary, we have analyzed a scheme for manip- 
ulating and measuring quantum states of single atoms 
in optical lattices with the assistance of focused lasers. 
With properly chosen experimental parameters, various 
manipulation errors are suppressed below 10 -4 except 
the spontaneous emission probability of atoms in focused 
lasers. Our proposal includes a realistic and practical 
quantitative analysis suggesting plausible experimental 
implementation of single atom manipulation in optical 
lattice quantum computation architecture. The tech- 
nique also has broad applications for investigating various 
interesting physics in optical lattices such as engineer- 
ing and probing many-body quantum states in strongly- 
correlated systems, topological quantum computation, 
etc. [2l]]. We believe that our work establishes a prac- 
tical and workable method for single atom manipulation 
in optical lattices. 
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